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Abstract 

The  chi-process  is  defined  as  a  natural  generalization  of  the   chi- 
distribution  of  statistical  theory,     A  formula  is  derived  for  the  expected 
number  of  level-crossings  per  second  of  the  chi-process.     The  formula 
contains  as  a  special   case  the  familiar  expression  for  the  fluctuation 
rate  of  the  envelope  of  Gaussian  noise. 
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1,  Introduction 

The  expected  number  of  times  per  second  that  a  stationary  random 
process  ^(t)  assumes  the  value  a  is  given  by  the  familiar  expression*-  -J 

00 


(1)       N(a)  -  J  \i\   p(a,^)  di, 

-00 


where  p(a,4)  is  the  bivcriate  probability  density  of  K  and  its  time  derivative  ^, 
evaluated  at  5  =  a.  Consider  first  the  case  where  C  denotes  a  zeixs-mean  Gaussian 
process  g(t)  with  the  correlation  function 


00 

R(ir)  =    <rg(t)g(t^r)>      ^     r  ^(^)^os  cordco. 

<g^(t)>         i 


Here  the  angular  parentheses  denote  the  expectation  or  ensemble  average,  and 
^(co)  is  the  (normalized)  power  spectrum  of  the  process  g(t).  Then,  equation  (1) 
reduces  to"-  -^ 


(2)        N  (a)  =  /I  (r/:^  W  (a), 


where  cr    =    <g  (t)>  ,  R"  is  the  second  derivative  of  R(1;)  with  respect  to  V 
evaluated  at  ^=  0  ,   and 

W  (x)  -  (1/ /2w  0-)  exp(-x^/2o-^) 

is  the  univariate  probability  density  of  the  process  g(t),  i.e.,  the  probability 

2 

density  of  a  Gaussian  random  variable  with  mean  zero  and  variance  <T    . 


Alternatively,  (2)  can  be  written  as 

'— \l/2 

where 


(3)      Ng(a)  =/|(y(co2)    Wg(a), 


*The  prime  will  be  used  to  denote  differentiation  with  respect  to  t ,   and  the 
dot  to  denote  differentiation  with  respect  to  t. 


-  2  - 


00 

T 

CO 


=  j   (0  ^(co)dw 


is  the  mean  square  (angular)  frequency  of  the  process  g(t).  Consider  next  the 
case  where  ^  denotes  the  envelope  E(t)  of  the  Gaussian  process  g(t).  Then, 
as  shown  by  several  authors'-  i'l^i'^J,   equation  (l)  reduces  to 

(k)  Ng(a)  =/|cr^coW2(a), 

where 

WgCx)  »  (x/a^)  e3q3(-x^/2(3-2),  x  >  0, 

(5) 

W-Cx)  =0  ,  X  <  0 

is  the  univariate  probability  density  of  E(t),  i.e.,  the  familiar  Rayleigh 
distribution.  The  quantity  ^co  is  given  by 

(6)      (^co)2  =  7    -  (SJ)2  , 

where  co    is  the  mean  square  frequency  of  the  underlying  Gaussian  process  g(t), 

00 


j       o>  J(c 


CO    »     I      o>  $(co)cia 
o 


is  its  mean  frequency.  It  follows  at  once  that  any  two  of  the  quantities 

S,  W  (a),  N_.(a)  determine  the  other.   The  similarity  of  (3)  aid  (U)  is  apparent. 

In  this  note  we  define  an  infinite  class  of  random  processes,  the  so- 
called  chi-processes.  We  find  that  the  expression  for  the  fluctuation  rate 
of  these  processes  is  a  natural  generalization  of  (2),  (3)>  and  (U). 

If  we  regard  J(co)  as  a  distribution  of  unit  mass  on  the  interval  (0,oo), 
then  a  is  the  centroid  of  ^(co),  ca  is  the  moment  of  inertia  of  J(co)  about  the 
origin,  and  (^<o)  is  the  moment  of  inertia  of  J(co)  about  its  centroid.  With 
this  interpretation,  the  relation  between  co,  N  (a),  and  N_(a)  is  an  instance  of 
the  parallel  axis  theorem  of  mechanics. 
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2,  The  chi-process  and  its  fluctuation  rate 

We  begin  by  defining  the  chi-process  by  generalizing  to  random  processes 
the  familiar  definition  of  the  chi-distribution"-  ■* . 

Definition;  Let  g-(t),,..,g  (t)  be  n  independent  stationaiy 
Gaussian  random  processes  such  that 

<  gi(t)>  -  0  , 
(7)  <§?(*)>  =  ""^  ' 

<  gi(t)gi(t+r)>  =(T^R(r),   1  <  i  <  n. 

Then  by  the  chi-process  ( 9^ -process)  with  n  degrees  of  freedom  is 


(8) 


meant  the  random  process 

L   1=1      -i 


The  processes  &,(t),  l<i<n,  can  be  regarded  as  the  components  of  a 
vector  random  process  g(t),  and  the  process- A  (t)  as  the  length  of  g(t).  The 
conditions  (?)  mean  that  the  components  of  g(t)  have  mean  zero  and  are  identi- 
cally autocorrelated.  The  univariate  distribution  of  the  process  7^^^*^  ^^*   °^ 
course,  the  X -distribution,  with  probability  density  W^(x)  given  by"-  J 


2    cr  I  t.-^; 


n^ 


(9) 


W  (x)  =  0  ,       X  <  0. 

n 


(Note  the  consistency  of  ($)  and  (9).)  We  now  derive  an  expression  for  the 
fluctuation  rate  of  the  chi-process. 

Theorem:  The  expected  number  of  times  per  second  that  the  process 

%   (t)  assumes  the  value  a  is  given  bj 


-li- 
do) N  (a)  =  /I    o-yCFi    W^(a). 
n  /    Ti        'on 

Proof.     We  consider  first  the  trivial  but  exceptional  case  n  =  1,   and  note  that 
Xi('t)  =  /g^(t)     =   |g(t)|.     Since   |g(t)l   »  a  >  0  if  and  only  if  g(t)  =  a  or 
g(t)  =  -a,   it  follows  from  (2)  that 

N;^(a)=/|a    /I^OTgCa),  a>0     , 

N^(a)   =0  ,  a  <  0     . 

Since  by  (9) 

W^(a)   =  2Wg(a)     ,  a  >  0, 

W^(a)   =0  ,  a  <  0, 

equation  (10;  is  established  for  the  case  n=l. 

When  n  >  2,  we  must  apply  the  basic  formula  (l),  which  requires  that  we 
first  calctilate  the  bivariate  probability  density  p{'%   >7^„)  o^  r-     a^ri  its 
derivative  %  •  We  begin  by  writing  the  probability  density  of  the  2n  jointly 
Gaussian  random  variables  x^«=g-, , .  ..,x  =g^,  x  ^-=g^,...,  Xp  =g^.  (Since  all  the 
processes  are  stationary,  the  argument  t  can  be  dropped.)  It  follows  from  (?) 
and  the  independence  of  the  processes  g.(t),  1  <  i  <  n,  that  the  matrix  M  with 
elements  M.  .  =  ^x.x.^  is  diagonal  with  non-zero  elements 

^-  •••  -^nn-"^^   > 

Vl,n.l"  •••=  V2n  =  -^'Ro  • 
The  probability  density  of  the  2n  random  variables  x,,...,X2     is 

"^Unlike  the  case  n  >  2,'/^(t)   assumes  the  value  zero  (without  passing 
through  it)  at  a  non-zero  rate.     y\lthough  the  sample  functions  of  ^(t)  are 
not  different iable  at  their  zeros,   the  latter  are  a  set  of  measure  zero,   so 
that  the  sanple  functions  are  differentiable  almost  everywliere. 
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where  M~  are  the  elements  of  the  reciprocal  matrix  M  .  Since 


det  M  =  ^^"(-R'y  , 
o       ' 


,-1 


and  M       is  diagonal  with  non-zero  elements 


•?i 


M 


-1  ^    1 
nn  "^  "12 


M 


.-1 


M, 


,-1 


n+l,n+l       •••       ''2n,2n  2l^     ' 

O 

we  find,   returning  to  the  variables  g.   and  g.,  that 


r     n 


P(g-j^»  •  •  •  yg^jg-j^j  •  •  .  *gj^) 


(2.<r2)"(-R^)"/2 


exp 


^     ^i  Z,     % 

1=1  1=1 


20-' 


5— rr 

-20-^  R^ 
o 


We  now  introduce  n-dimensional  spherical  coordinates  by  writing 


%       "  9^n  ^^"  ®1  ^^'^  ^2  • "  ^^"  ®n-l  * 
gg       "   fn  ^^  ^1  ^-^^  ©2  •  •  *   °°^  ®n-l 


Sn-l  -  /n  ^^  ®1  °°^  ®2 


^n       -/n 


cos  9, 


which  maps  the  region  -oo<g.  <oo,  liisn,  into  the  region 
^  -f-n^  00, 
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0  ^  ®i»*'»'®n-2  -  "» 


°  ^  ^n-1  <  2n. 


We  see  at  once  that 


and  after  a  bit  of  manipulation  that 
n 


(11)    Z  ^?    -X?  -7^^  *?*  f  ^  i\  3in\  *  ...»  ^^  9^1  -A  •■ 

1=1 


.  sin' 


h 


n-2' 


Thus,  the  joint  probability  density  of  the  new  random  variables  ^  ,9^,..,,©       , 


P^'^n'®l»  •••»®n-l'?^n'^l'  *  *  "^n-l^ 

i^  n  ?^  n*^n^l*7-^2^i"  V '  *  -^Cl^^^^- ' ' ^^"^-2 


(2n<r2)^-R")°/2 


exp 


2a' 


-20^  r" 
o 


S(g]^»g2>»«»>gj^>gl>g2»****^n^ 

*     *  * 

S(  "X-j^J^i*  •  •  •  >®n-l»  ^n'^1'  •  •  •  »^n-l^ 


,W 


It  follows  from  Laplace's  expansion  theorem  for  determinants*-  J,  and  the  fact 

that  the  partial  derivatives  of  the  variables  g-i,g«,...,g  with  respect  to  the 

variables  7-  j®n*»«»i9  t  vanish,  that 
/  n  1    '  n-i      ' 


/  *  *     *  \ 


9(g2^>g2*«»«*gn)    S(g^,g2,...,g^) 


3(  ;^n'®l*  •  •  •  '®n-l»  ^n*®l*  *  * '  '^n-1^    ^^^  n*®l»  *  * '  '®n-l^  ^^  "^n'^l*  *  *  *  »^n-l^ 


^iz-^i%2>  *  *  *  »^r) 


e(Xn'®i,...,Vi^ 


n 
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where  the  Jacobian 


J  =  l^n-l  sin^^-^e  sin"~%-  ...  sin  9  - 
n    n        1       <:        n.*^ 


is  the  square  root  of  the  product  of  the  coefficients  of  the  quantities 

OC^,§T,.,.,e^  T  in  the  svun  (11). 
n'  1'   '  n-1 

The  bivariate  probability  density  p(  'f-  y  %  )   is  given  by 
n     n    2n    00     CO 

p(T^^,  i^)   -  /dftL-./«r>.2Pn-lf^r--f  ®n.lP^^n'®l--^-l'  Vl--^-l^ 
0     0    0    -00    -co 


which  equals 

n 

(12) 

00 

* 
X    c3&,exp 
J 

-00 

\   i-Yr  1 

dOgexp 

^00 

-2a'^R 
o. 

exp 


2     •  2   -[ 

^n     "^  n 

— 5 5~n- 

2ct'^   -2a'^R  ^ 
o  -* 


2a2 


.iZi 


sin 


■^. 


-2cV 


CO 

...j  d9^_^exp 

-00 


"2*2     2       2« 
-^  9  ,  sin  ©, . . .  sin^  « 
Xn  n-1    1 n-2 

-2(r^R" 

0 


Tl^ie  integrals  with  respect  to  ®i>»*'»^n-l  *^^  ^®  ^°"®  immediately.  Their  product 


is 


2  " 

o 


% 


n 


n-1 
"5" 


sin""  9,sin'^"-^9^...sin  9„  « 
1      2       n-£ 


n-1 
(.2.o-2lfe")"^ 


n 


so  that  (12)  becomes 


P<^n'  T'n' 


1 


(2n<r^)^^  /^^     o  o 


n  It  2n 

[d9^...|de^_2(d&^_^  J„exp 


^2      ^2     n 

^n       -^1 

2a^    -2cT^"   . 
o  -• 


Since 


n  n  2n 
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where 


2n 


n/2 


r(f) 

is  the  area  of  the  n-dimensional  sphere  with  unit  radius,  we  find,  using  (9),   that 


^  r n'  'n  n   '  n 


'  0 


exp 


'    2 

'^  n 
5-n- 

I.  o     J 


It  follows  that  at  any  time  t,  the  process  '^  (t)  and  its  derivative  'X  (t)  are 


n 


n 


independent.  Indeed,  the  univariate  distribution  of  'X  (t)  is  Gaussian,  with 


2  "  w  /  \ 

mean  zero  and  variance  -o"  R  ,  although  its  integral  'J-   (t)  is  manifestly  non- 
Gaussian,  This  involves  no  contradiction,  since  the  process  7^  (t)  is  not 
multivariate  Gaussian.  The  proof  of  the  theorem  is  completed  by  noting  that 


00 


n 


00 


exp 


-co 


•  2 

An 


2 — 


'A. 


o  -^ 


so  that 


(10)      Nj^(a)  =  /|  o/~R^\i^U),       Q.  E.  D. 

In  proving  the  theorem,  the  independence  of  the  processes  g-.  (t),,..,g^(t) 
was  used  only  to  infer  that 

(13)   <gj^(t)g^(t)>  =  <gj^(t)ij(t)>-  <q(t)g^(t)>  =0,  1  <  i  /  j  <  n. 

(The  relations  <Cg.  (t)g.(t))>  =0,  1  <  i  <  n,  follow  automatically  from  R'(0)  =  0.) 
Thus  we  have  the  following  result. 


Corollary;  Let  g-,(t),..,,g  (t)  be  n  jointly  Gaussian  processes  which 
are  stationary  and  stationarily  correlated  and  which  satisfy  (7). 
Suppose  in  addition  that  the  cross-correlations  of  g_ (t),. .,,g^(t) 
are  such  that  the  conditions  (13)  are  satisfied.  Then  the  expected 
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rl/2 


X  g?(t) 


assumes 


number  of  times  per  second  that  the  process 

the  value  a  is  identical  with  that  of  the  '^  "-process  ^  ai^d  is  given 

b2  (10). 

3,  Application  to  the  envelope  fluctuation  rate 

We  now  show  that  the  general  expression  (10)  includes  as  a  special  case 
the  formula  (h)   for  the  fluctuation  rate  of  the  envelope  E(t)  of  a  zero-mean 
stationary  Gaussian  process  g(t).  First  we  recall  the  following  facts  about  the 
envelope  1-  -• . 

1)  For  every  co  >  0  there  is  a  representation  of  g(t)  of  the  form 

g(t)  =  g^(t)  cos  co^t  -  gg(t)  sin  co^t, 
where  the  processes  g.(t)  and  g^(t)  are  jointly  Gaussian,  and 


<g^t)>     -    <g3(t)>    =    <r2. 


(lU) 


where 


(15) 


2)  The  correlation  properties  of  g  (t)   and  g„(t)  are  given  by 

C  3 

<gc(t)gjt+r)>    -    <gg(t)gg(t*t-)>    -cr^RCr), 
<  g^.(t)gg(t*  r)>    = -<g^(t)g^(t+r)>  -  <j2s(t-). 


00 


R(T)  =   I    _|(co)  cos[(co-«)^)  r)d£o. 


o 

CD 


S(r)  =  j  J(aj)  sin[(co-oj^)  rjcko. 
o 
The  quantity 

(16)  k2(r)    -   R^C^i  +  s^(i-) 
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is  independent  of  co  , 


4* 

3)     In  particular,  it  follows  from  (Hi)  and   (15)   that 
<g,(t)g^(t)>     -     S^     -     0, 

<gc(t)g,(t))  -    S^   -    0, 

and  that 


doi     "  CO  -  0)     , 
o   ' 


where  w  is  the  centroid  of  J(oo),  regarded  as  a  mass  distribution. 

h)     The  envelope  E(t)  is  defined  by 

1 

E(t)   .   [g^(t)  .   g^(t)J      , 

and  is  independent  of  the  choice  of  w  . 

It  follows  that  the  conditions  needed  to  apply  the  corollary  of  the 
fluctuation  rate  theorem  are  met,  provided  that  we  choose  the  frequency  oo  , 
which  is  at  our  disposal,  to  be  the  centroid  of  the  distribution  J(co).  With 
this  choice 

-00 


-  R^'    -  (oi-«)    J(w)  dw     "     (Acu)^   , 


v2 


where  (Aoo)     is  the  quantity  defined  by  (6).     Thus,   the  envelope  fluctuation 
rate  N„(a)  is  given  by 


N^(a)     -     NjCa)     -    /f  '^/^  ^2^®^' 


or 


The  subscript  eero  denotes  evaluation  at      6    ■     0. 
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(U)  Np(a)  -  y/T'<rAcoW2(a), 

which  completes  the  proof  that  (U)  is  a  special  case  of  (10), 

It  is  interesting  to  note  another  expression  for  if-c.{a)   in  terms  of  the 
envelope  correlation  function 

p  (r)   .   <E(t)E(t.r)>  -  <E>^      _ 


It  can  be  shown  that 


[7] 


R  (r)   -  ^^(^)"(i->t^)K(^^)-("A) 

^  2  -  (n/2) 

where  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind,  respectively,  and  k  is  the  fvinction  of  t'  defined  by  (16). 
It  follows  that 

(17)     -  R^(0)  -  jl^   0/3k)J2E(k)-(l-k^)K(k)l   (-R^"  -  s'^)  -  jL.   (Ao))^  , 

k-1 

where  we  have  used  formulas  112,01,  710.02,  and  710 ,0U  from  Byrd  and  Friedman's 
handbook  of  elliptic  integrals'-  -I ,  Comparing  (U)  and  (17) i  we  see  that  another 
expression  for  N_(a)  is 

Ng(a)  -  ^/^  <r^-R^'(0)  V^{b). 
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